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V/-^ ■ Abstract: It is very interesting that all holographic superconductors, such as s-wave, 

^ ■ p-wave and (i-wave holographic superconductors, show the universal mean-field critical 

exponent 1/2 at the critical temperature, just like Gindzburg-Landau (G-L) theory for 
second order phase transitions. Now it is believed that the universal critical exponents 
appear because the dual gravity theory is classic in the large limit. However, even in 
^ ■ the large N limit there is an exception called "non-mean-field theory" : an extension of 

■ the s-wave model with a cubic term of the charged scalar field shows a different critical 

exponent 1. In this paper, we try to use analytical methods to obtain the critical 
exponents for these models to see how the properties of the gravity action decides 
the appearance of the mean- field behaviors. It will be seen that just like the G-L 
theory, it is the fundamental symmetries rather than the detailed parameters of the 
bulk theory that lead to the universal properties of the holographic superconducting 
phase transition. The feasibility of the called "non-mean-field theory" is also discussed. 
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1. Introduction 

Using AdS/CFT correspondence |1|, ^ |], Q to study superconductors is a very inno- 
vative idea. The method to build the holographic superconductors that we use here 
is called the "bottom-up" approach: by putting some set of fields in the asymptotic 
AdS spacetime, through the study of the weak interacting bulk field theory we can get 
the information of the strongly coupled boundary theory. For example, when putting 
a charged scalar field coupled to a f/(l) gauge field in the bulk, we get the s-wave 
holographic superconductor in which the order parameter is a scalar]^, ^. By putting 
a pure SU (2) gauge field in the bulk, we can get the p-wave holographic superconduc- 
tor with a vector order parameter |0, ^, ^. In order to build a d-wave holographic 
superconductor we need a charged tensor field coupled to a gauge field in the bulk 



that leads to a tensor order parameter |]10|, 0, |12[- Besides the "bottom- up" approach, 
there is also a "top-down" way to embed the holographic superconductors into string 
theory, since the scalar, gauge and tensor fields we used above have a natural source 
in string theory. In Ref. the authors studied various D-brane configurations for 
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the p + ip-wave and s-wave model, whereas the embedding of the d-wave model with 
tensor field into string theory is still an open question . 



These holographic superconductors have many properties similar to real world su- 
perconductors: the second order phase transition at the critical temperature, the be- 
haviors of holographic superconductor under a magnetic field, the AC conductivity, and 
so on. One of the most important success of holographic superconductors is that all 
these models show the exact mean-field behaviors at the critical temperature, just like 
the Gindzburg-Landau (G-L) theory for continuous phase transitions. All the critical 
exponents for the order parameter at are 1/2. The G-L theory is a tremendous 
breakthrough towards understanding the physical spirit of the continuous phase transi- 
tion which enables us to understand that it is the fundamental symmetry that decides 
the universality of continuous phase transitions. Till now, the study of holographic 
superconductors made us believe that the AdS/CFT correspondence indeed contain 
the physics of real world superconductors. We know that the fundamental symmetry 
decides the rule of the G-L theory. It is also a meaningful thing to answer whether this 
also works in the holographic models. However, most of the computations to obtain 
the critical exponents are numerical. It is hard to see exactly what properties of the 
action decide the mean-filed behaviors. In this paper, we employ the analytical method 



14| to stuty the s-wave, p-wave, d-wave and extended s-wave holographic superconduc- 



tors. There are other papers that study holographic superconductors with analytical 



methods p5|, P4| . The analytical calculation indeed provides more information to an- 
swer where the universal critical exponents 1/2 of the order parameter comes from. 
It is indeed the fundamental symmetry of the bulk theory that leads to the expected 
mean-field behaviors. 

The organization of this paper is as follows. In Section 2 we analytically compute 
the critical exponents at for s-wave, ci-wave and p-wave models, respectively. We 
can see how the fundamental symmetry of the actions decide the critical exponents. 
In Section 3 we give the computation of the extended s-wave model with mean-field 
behavior and give some comments on the critical exponents different from 1/2 in the 
extended s-wave model. Discussions and conclusions are given in Section 4. 



2. The s-wave model, j9-wave model and d-wave model 

The construction of holographic superconductors is based on the fact that a black 
hole coupled with charged matter fields may have hairs and then will break the local 
U(l) gauge symmetry in the bulk. This happens when there is an instability towards 
a charged matter field to get an expectation value due to the bulk action and its 
boundary condition. To construct different kinds of holographic superconductor, we 
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use different matter fields in the bulk action and solve the equation of motion to get the 
background. Generally speaking, this can be done only numerically. However, there 
always exists some parameter space which allows us to decouple the matter fields and 
the AdS black hole background. This decouple limit which we will use in the following 
sections is always taken place in the large charge limit of the matter fields and is called 
"probe limit" . Firstly, we introduce the AdS black hole background that will be used 
for several holographic superconductor models. The AdS black hole background can 
be used only if it is a solution of the equations of motion. Then we will analytically 
compute the critical exponents for these models. 

2.1 The Anti-de Sitter black hole background 

According to the AdS/CFT dictionary, if we want to study a finite temperature field 
theory on the boundary, we need an AdS black hole background. The added matter 
fields are perturbation of the black hole in the probe limit and the metric will not be 
affected by the added matter fields. The metric of the 74^5*4 black hole in units in which 



the AdS radius L = 1 reads 14 



ds^ = -f{r)de + r^dx" + dy^) + (2.1) 

fin 



where 



2 



f{r)=r'-^. (2.2) 
r 



The Hawking temperature of the black hole is given by 

T=^±, (2.3) 

which is also the temperature of the dual gauge theory living on the boundary of the 
AdS spacetime. This is a AdS^/CFT^, correspondence, which means that the dual 
superconductor is of 2+1 dimension. In the probe limit, holographic superconductors 
are constructed by adding perturbative matter fields to the black hole background. For 
example, for the p-wave model we need SU{2) gauge fields in the bulk. 

2.2 the s-wave model 

In this section we review the first analytic computation of the critical exponent of the s- 
wave model which has been made in Other analytical methods have been developed 



in[pl|, p4| . For review of the s-wave model one can see |1l6| , 0, |T8|, 0. The matter 



fields of the s-wave model include a charged scalar field and a Maxwell field coupled to 
it. The condensation of the scalar operator on the boundary which corresponds to the 
scalar field in the bulk is the order parameter of the superconducting phase. 
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= -^F^^F^"" - - iA^p - m^l^l^ (2.4) 

We take the ansatz that ^(r) 7^ and $(r) 7^ 0, $(r) is the scalar potential 0. 
We introduce a new coordinate 



r 



(2.5) 



The equations of motion (EOMs) in the coordinate z are 



3 



2;^ - T r^' + 

1 — Z'^ 



$2 ^2 



'r2 (1-^3)2 ^(1_23) 



^ = (2.6) 



where the prime denotes differentiation with respect to z. We have to solve the EOMs 
in the interval (0, 1), where 2; = 1 is the horizon and z = is the boundary. As will 
be seen below, by solving the equation of \E', one can obtain the value of Tc, while by 
solving the equation for $, one can find the critical behavior of the order parameter at 
Tc. Near the boundary (2; — )■ 0), the boundary conditions are 



^-^-^^ , $-/.-^. (2.8) 



where 



A± = ^±^^ + m2 (2.9) 

(Oa±) are the condensation with dimension A-t, is the chemical potential and p is 
the charge density of the field theory. 

At the horizon, we require that the scalar field be regular and 

<I'(1) = (2.10) 

to ensure At = ^dt to be well defined. 

Since there is a second order continuous phase transition at the critical temperature, 
the solution of the EOMs at Tc is 

^z) = \r+,{l-z) , A = -^ (2.11) 
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where r+c is the radius of the horizon at T = Tc- 

As T — )• Tc, the field equation for the scalar field \l/ approaches the limit 



2 + 



m 



1 + z + z' 



2\2 



(2.12) 



By solving the equation of \E', we can obtain the value of T^. To match the behavior 
at the boundary, we can define 



z'^F z) 



(2.13) 



where, according to eq. ( p.8| ), F is normalized as F{0) = 1. 
We deduce 



-F" + 



2 + z^ 
1^ 



-2A 



m2 + 3A - A2 + A2^3 



(1 + z + 2;2)2 



F (2.14) 



when z — )■ 0,^ should be finite, so this equation is to be solved subject to the boundary 
condition 

F'(0) = (2.15) 

In order to compare the analytical results to the numerical results in P, we set 
= —2, then from eq.( |2.9| ) the dimension A of the order parameter (Oa) are A = 1 
or A = 2. By putting = —2, the equation for F reads 



F" + - 

z 



2 + z^ 



2A 



F' 



A^z 



l + z + z' 



2\2' 



(2.16) 



Now, eq.( ^.16 ) is translated to the Sturm-Liouwille form under the certain bound- 
ary condition. According to the general variation method to solve the Sturm-Liouwille 
problem (|^T|,or see Appendix), the eigenvalue minimizes the expression 



A^ 



/g dz ^2A-2|(i _ z^)[F'{z)]^ + A^z[F{z)]^} 



(2.17) 



In order to use the variation method , we have to specifies the trial eigenfunction F{z). 
From the boundary condition F{0) = 1 and F'{0) = 0, the constant term is 1 and the 
linear term of z is forbidden. The second order trial eigenfunction is then 



az 



(2.18) 
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If the lowest solution of F{z) corresponds to a > 0, the whole solution of ^{z) will 
decrease monotonically from the boundary to the horizon and it is believed to describe 
the superconductor vacuum solution |]19[. It will be seen in the following that this is 
true. 

For A = 1 we obtain 

= j ' + (2.19) 

2y37r - 6 In 3 + 4( VStt + 3 In 3 - 9)a + (12 In 3 - 13)a2 

which attains its minimum at a ~ 0.239. We obtain 

A' ~ ^ 1.268 (2.20) 

which can be compared with the numerical value A^ = 1.245 0. The critical tempera- 
ture is 

Tc = — r+, = — ./^ (2.21) 

so for A = 1, Tc ~ 0.2250y^, which is in very good agreement with the numerical value 
T, = 0.226yp i. 

In fact, this analytical calculation can be done even better if we include higher 
order of z such as third order trial eigenfunction F = F^islz) = 1 — az"^ + Pz^. Then 
the eigenvalue A^ minimize eq.( p.l7|) when a — ?• 0.4962, /3 — )■ 0.2706 . As a result 



A^ ~ 1.259 and Tc ~ 0.2255^//) , which is more close to the numerical result. However, 
for qualitative analyze, the second order trial eigenfunction is good enough and we use 
it in the following sections. 

Similarly, for A = 2 we obtain 

A^ = 2 l-la + ja^ ^ 

3-ln3-^ + (f -41n3)a + (^- j^ + ln3)a2 v ; 

whose minimum is A^ ^ 17.3 (at a ~ 0.6), which can be compared with the exact value 
A^ = 16.754. The critical temperature in this case is ~ 0.117 ^/p, which is in very 
good agreement with the numerical value Tc = 0.1 18 0. 

Now we begin to solve the equation for $ to obtain the behavior of the order 
parameter at Tc. Away from (but close to) the critical temperature, the field eq.( p.7|) 
for $ becomes 

(Oa)^ z'^^'^-'^'>F^(z) 
*" = ^ ' * (2-23) 
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where the parameter {O ^)'^ / (2r']^) is small. We may expand $ in the small parameter 



as 



- = A(1-^) 



(2.24) 



where x is the general correction function. For the correction x near the critical tem- 
perature, from eg. ( p.23| )we deduce 



X" = A- 

with x(l) = = 0. 

To find the temperature, we need 

X'(0) = AC , C 

From eg. ( |2.8| ), we deduce the ratio 



4 = A(1 



(2.25) 



dz 



1 + ^ + ^2 



(2.26) 



^2A 



+ . . . 



(2.27) 



Using the fact that = "i^J p/X/A-n and T = 3r+/47r, the eguation becomes 



Then 



T^^-^(r, + r)r,(i--)oc(OAr- 



Since T is very close to Tc, we have 



T 



1/2 



(2.28) 

(2.29) 

(2.30) 



From the computations from eg.( |2.23| ) to eg.( |2.30| ) we can see that the critical 



exponent 1/2 comes from two aspects: the first one is oc a/ p/A, the second one 
is the \E'^$ term in eg.( |2.23| ) which determines the order of the order parameter in 



eg. ( p.27| ) to be 2. The fact that Tc oc \/~pJ\ is from dimensional analysis, while the 
second order of (0a) in eg. ( p.23|) is determined by the ^1^^$^ term in £s, and this term 
is determined by the fundamental symmetry. In Ref. [p5[] the author also found the 
results that the critical exponent 1/2 follows from the \E'^$ coupling. The fundamental 
symmetry includes the local t/(l) gauge symmetry and the positive/negative symmetry 
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of the order parameter. The latter one determines that there are only quadratic and 
quartic terms in the action. As we will see in section 3, the critical exponent differs 
from 1/2 when the added cubic term breaks the positive /negative symmetry of the order 
parameter. If one adds quartic terms in the action, the value of the critical exponent 
will not be altered. The cubic term is not allowed in the G-L theory, since we have 
to preserve the positive/negative symmetry of the order parameter as a fundamental 
symmetry. If we add the cubic term in the G-L theory anyway, we will naturally get a 
critical exponent different from 1/2. If we put the ansatz \E'(r) 7^ and $(r) 7^ into 
the action and do the computation, we will get 

= ^ fi^d^mf + ^2^2 _ ^Q^^Y ^2.31) 

This Lagrangian has the same symmetry as the G-L Lagrangian. 
2.3 the p-wave model 

The matter field of the p-wave holographic superconductor is a pure SU(2) gauge field 



= -^TT{F,^Fn. (2.32) 

where qym is the gauge coupling constant and F^y = T°'F^^ = d^j_Ay — dyA^ — i\A^, A^] is 
the field strength of the gauge field A = A^dx^^ = T"-A'^^dx^. For the SU (2) gauge sym- 
metry, [T",T^] = ie°'^'^T^ and Tr (T^T'') = 5"''/2, where e"'^'^ is the totally antisymmetric 
tensor with e^^^ = 1. The Yang-Mills Lagrangian becomes Ti{F^yF^^) = F^^F"-'^^ /2 
with the field strength components F'^^ = d^A^ — d^A'^^ + e^-'^^A^^A^. 

The ansatz for the p-wave background is A = (f){z)T^dt + w{z)T^dx. with other 
components vanishing. Here the U{1) subgroup of SU{2) generated by is identified 
to be the electromagnetic gauge group and is the electrostatic potential, which 
must vanish at the horizon for the gauge field A to be well-defined, but need not vanish 
at infinity. Thus the black hole can carry charge through the condensate w, which 
spontaneously breaks the U{1) gauge symmetry. The EOMs for the two fields (p and w 
are 



. ^_-,^ = (2-34) 

It is clear that the coupled term of the condensed matter field and the scalar 
potential is in which the condensed matter field w is quadratic. This is similar to 
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the s-wave model which is important to obtain the critical exponent 1/2, as we have 
discussed in the last section. Now our task is to solve the EOMs in the interval (0,1) 
under the proper boundary conditions. The boundary condition at 2; = is 

w!=^^z, (1)^11-— z. (2.35) 

O is the order parameter, fi is the chemical potential and p is the charge density of the 
field theory. In the s-wave and rf-wave model, the dimension of the order parameter is 
determined by the value of the mass. While in this p-wave model the order parameter 
is just of dimension one. At the horizon, we require w be regular and = 0, since (j)dt 
needs to be well defined at 2; = 1. 

If there is a second order continuous phase transition at the critical temperature, 
the solution of the EOMs at the Tc should be 

w{z) = 0,cj){z) = Xr+,{l-z), (2.36) 

where A = -§-, r+c is the radius of horizon corresponding to T" = Tc. So, at a temper- 
ature slightly below Tc, the EOM for w becomes: 

z{i — z'^) [1 + z + z^y 

We can define w as 

w(z) = ^zF(z) (2.38) 
r+ 

To math the boundary condition at 2; = 0, i^(0) should be normalized as i^(0) = 1. 
Then we can get the equation for F{z), 

,, 1/ 3zF(z) ^ ^ ^ 

-F" + -( --2)F' + - = 7 2.39 

zh-z^ ' {l-zf (1 + ^ + ^2)2 ^ ' 

By solving this equation, we can obtain the value of the critical temperature. Besides 
the boundary condition F(0) = 1, another boundary condition is -F'(O) = 0. Then we 
can solve this equation by using the same method for the s-wave model. The eigenvalue 
A minimizes the expression 

^2 ^ Ipdz zmi-z')[F'{zr + 3z[F{z)n 

To estimate it, we use the trial function 

F = F^{z) = \-OLZ^ (2.41) 
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We then obtain 

\^ 10 ('2 42) 

" 6-21n3-^7r-(81n3-f)a-(| + 21n3-^7r)a2 

which attains its minimum at a ~ 0.51. We obtain 

^ 51 ^ 13.77 (2.43) 

The critical temperature is 

So, Tc ~ 0.124y^, which is in agreement with the numerical result Tc ~ 0.125^/p in 
i- 

If we want to know the behavior of the order parameter at Tc, we need to solve the 
equation for the scalar potential close to 

^ = ^2 1 _ ^3 (2-45) 

Since the order parameter {O) is small, is a small parameter. We can expand in 
this small parameter as 

- = A(l-^) + ^xW + ---. (2.46) 

r+ ri 



Then we get the equation for x. 



X" = A ' V ^2.47) 

1 + 2; + 



and x(l) = X'(l) = 0. 

From the asymptotic behavior of (j) at the boundary, we have </''(0) = p/r^. We 
can also get x'(0) from the above equation, 

X'{0) = XC,C= dz \\ . (2.48) 
Jo i + z -\- z 

We have 

'+ ' + 

Then 

(O) oc Tc(l - ^y/', (2.50) 
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which is the exact mean field critical exponents given by numerical calculations. The 
]9-wave model is different from the s-wave and ci-wave model. In this p-wave model the 
condensed charged matter field is not coupled to the U{1) via the covariant derivative. 
It is hard to reveal the fundamental symmetry of the p-wave action. If one puts the 
ansatz into the action, it is clear that the action has the same fundamental symmetry 
as the s-wave model with terms like (dzw)"^, (5^0)^ and w'^cf)^. 

2.4 the (i-wave model 

In order to construct a holographic model of d-wave superconductor, we need a spin 
two field in the bulk [0. In JTT], 0, a d-wave model with charged tensor field and 
Maxwell field in the bulk is built, in which the inconsistency of charged tensor field in a 
curved background like the appearance of ghost and causality can be made very small 
in some limit. The matter field of the holographic d-wave model is. 



(2.51) 



where = — iqA^ and ipp = D^ip^p. The equations of motion which follow from 



are 



= (□ - m^)ipp^ - 2D(^pipy) + D(^pD^)ip - g^y [(□ - m^)ip - D^'ipp] 

+ 2Rppy,ipP^ - 9,u^^ - i\ {F,p^t + Fyp^^^) (2.52) 

D^F^^ = r 

where 

r = tipl^iD^if'^^ - D'^cp"^) + liifl - D^ip*){ip'"' - g'^^cp) + h.c. . (2.53) 
For the ci-wave backgrounds, the ansatz takes the following form |TT|, O 



l2 

A = Ap dx^ = <f){z) dt , ij,y{z) = — ^{z) , (2.54) 

2^^ 



with all other components of ipp^y set to zero, and (p and ip are real. The ansatz ( p.54| ) 
satisfies ip = ip^ = Fp^pipP = 0. Instead of turning on i/j^y in ( p.54| ) we could consider a 
non- vanishing value for ipxx-yy = i^xx = —i^yy These two ansatzs are equivalent under 



a 7r/4 rotation |T2 



With this ansatz, we can derive the equations of motion for the two fields ip and 
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2 -I- 



m 



^ = 



(2.55) 



It is interesting that the EOMs for the condensed field and scalar potential are the 
same as the ones for the s-wave model except for the coefficient of the term tp'^cj) in the 
equation for (f). But the boundary conditions at the boundary are slightly different, 



where p is the charge density of the boundary field theory, A is given by m'^L'^ = 
A (A — 3), and > ||12|. We can check that such a choice of dimension is compatible 



with the unitarity bounds in conformal theories, for instance In addition, in Ref. 
r^l, the authors give a generalization of the analysis of Breitenlohner and Freedman 
2T| to obtain the bound of > 0. Since the EOMs of the d-wave model are the same 

to those of the s-wave model and the boundary condition with different A will not 

affect our results of the critical exponents, we can conclude that for the d-wave model, 

the order parameter behaves as 

(a.) oc r,^(i - ^y/' (2.58) 



3. The extended s-wave model 



In order to construct a holographic superconductor with critical exponents different 
from 1/2, we need to extend the s-wave model by preserving the gauge symmetry. 
First we can rewrite the s-wave model innocuously in a Stiickelberg form by rewriting 
the charged scalar field as \E'e*^: 

= —F^uF"" - d^^ - ^^{dp - Af - m'l^l^ (3.1) 

So far all we have done is to rewrite the model. Nevertheless, it is straightforward to 
generalize the model in a gauge invariant way. The generalized action reads 

= -^F^uF'"' - d^^ - |/C(^)|(ap - Af - m^l^l^ (3.2) 
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This is called the Stiickelberg Lagrangian|2^, The general form of /C(^E') is 

/C(^) = ^C3|$|=' + C4|^|^ (3.3) 

In the above equation we have taken the absolute value since we require the Lagrangian 
still be local U(l) gauge invariant. In this model, the probe limit also works in the 
large charge limit ^3|]. When C3 and C4 vanish, the model reduces to the original s-wave 
model. When C3 = and < C4 < 1.4 |22|, the model has a second order phase 
transition with a critical exponent 1/2 at the critical temperature. Since this is still a 
continuous phase transition, we can apply the analytical analysis, as will be discussed 
in detail below. When C4 > 1.4, the superconducting phase transition turns to be of 
first order and our analytic method breaks down. 

With the ansatz $(2;) 7^ 0, 7^ and the gauge freedom to fix p = 0. we have the 
EOMs: 

2 + ;z3_ $2 



1 



+ Z 



2rl(l - z 



3^2 



m 

z{\ - z^ 
2/C 







■$ = 



Z\\ - Z3) 

The boundary condition for and $ are the same as that of the s-wave model. It is 
straightforward to repeat the computation for s-wave model here with also = —2. 
Now, the equation for F becomes 

"2 + ^3 



1-^3 



-2A 



1^ 



2^2 



(F + 2C4F'' 



(3.4) 



Since F is a small near the critical point, it is reasonable to ignore the 2c4F^ term 
if we are only concerned about the behavior near the critical point. After solving the 
equation of F with the same method used before, we obtain Tc ~ O.llTy/p for A = 2. 
For different values of C4, the critical temperature is the same. This conclusion is in 



agreement with the numerical computation in Ref. |22 



We now solve the equation for $ to see if the mean field behavior will be altered 
by /C. The equation for $ close to Tc is 



$"= ( 



^2A 



z^'^F^iz^ 



+ 



C4(Oa)^ 



z^^FHz\ 



2r 



4A 



(3.5) 



where C4 < 1.4. Similar to the computation in the above sections, we expand $ in the 
small parameter -^^j- as 



- = A(l-^) + ^x(^) + .... 
r 1 rl 



(3.6) 
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We can get the equation for x 

X = 2n / , 2. (^'^'^ + ^^z'^F')- (3.7) 
In order to match the asymptotic behavior of $ at the boundary, we have 



.2 -A(l + ^^ + ^;^ + ---), (3.8) 



where 



- Jo z\l + z + z^y - io ^ z\l + z + z^) ^^-^^ 
Solving eq.(p.8p for (Oa) and selecting the only physical solution with (Oa) > 0, we 



get 

(Oa) ocT,^(l-r/T,)i/l ^3_^0) 

This is in agreement with the numerical results in [p^ . 

The critical behavior with critical exponent 1 appears when the |\E'|'^ does not 
vanish. The situation with C3 = — 1 and C4 = 0.4 was studied numerically in ||22|, ^ with 
the critical exponent being 1 rather than 1/2. The order parameter has both positive 
and negative solutions since the cubic term breaks the symmetry between positive 
and negative see Fig. (2) in [^. Since a negative order parameter is unphysical, 
we must get rid of the solution with negative values. The analytical method we use 
here is based on the fact that when we lower the temperature from the normal state 
with the condensation vanishing to Tc, the order parameter goes continuously from 
zero to a finite value at the critical temperature. While for the situation with C3 = — 1, 
numerical calculations tell us that if we lower the temperature from a high value, we get 
a negative order parameter before we get the finite positive physical order parameter. 
So the analytical method we use here are not feasible now. However, just as was shown 



in Ref. another analytic method is developed to find the relationship between 



the value of critical exponents and /C. The cubic term breaks the positive /negative 
symmetry of the order parameter. This is not allowed in the G-L theory, since we have 
only squared and quartic terms which preserve the positive/negative symmetry of the 
order parameter. 

Since the mean-field behavior is protected by the fundamental symmetry of the 
G-L theory, even though different from the original s-wave model, the extended s-wave 
model contains a quartic term, we still have the same mean field critical behaviors. 
Therefore, if we add a cubic term in the s-wave model to break this symmetry, the 
appearance of so called "non-mean field" critical exponent which differs from the usual 
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value of 1/2 can be understandable. However, the model with non- vanish C3 has problem 
when we investigate the target space of the model. In (3.1) the action of a complex 
scalar is rewritten in terms of real fields \l/(its modulus) and p (its phase). In (3.2) the 
kinetic term is modified: the action is still gauge- invariant, but it describes a cr-model 
whose target space is parametrized by = (^I/,p) and it is in general not flat. The 
kinetic term of the a model is 

-g,,d^<fW<f^. (3.11) 
Using the expression (3.3) the metric is 

g,j = diag(l, ^2 + C3^=^ + C4^"). (3.12) 

When C3 = C4 = 0, it is the flrst holographic superconductor in 0. It is clear that 
the metric is flat, so the cr-model just describes a complex scalar fleld. For the case of 
C4 7^ 0, C3 = 0, the target space is no longer flat: its scalar curvature is 

(1 + C4^2)2 V ; 

The target space is curved, but it is still a smooth two-dimensional manifold with 
no singularities, so the model makes perfect sense. However, the story is different 
when turning on C3 (and setting C4 = 0, as turning on C4 as well does not change the 
conclusion) the scalar curvature is 

C3(4 + 3C3^) 

When \1/ — )■ the scalar curvature diverges. Therefore the target manifold has a 
singularity at the origin. In the classical theory, we could make sense of the a-model 
by removing the singular point. However in the holographic model this looks quite 
awkward, because above the critical temperature the classical solution is precisely \l/ = 
0, that is it sits at the singular point. This indicates that the theory with C3 7^ have 
problem.^ However, we worked in the probe limit, the theory may make sense when 
we include the back reaction of matter fleld on the metric of the AdS black hole. 



4. Discussion 

In this paper, under the large N limit when the quantum fluctuation is suppressed, 
we analytically study the mean fleld behaviors of the order parameter at Tc for four 

^The discussion of the feasibility of the extend s-wave model is based directly on the referee's 
argument. We thank the referee for pointing this important fact. 
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different holographic models of superconductor and it is found that these four models 
have similar properties. For each model there is a charged matter field coupled to a f/(l) 
background gauge field. After solving the two non-linear coupled EOMs with proper 
boundary conditions, the information we need for the strongly coupled boundary theory 
can be obtained from the asymptotic behaviors of the matter field and gauge field near 
the boundary. Numerically it is hard to see what properties of the bulk theory leads to 
the mean field or "non-mean-field" behaviors. However, by using the analytical method, 
it is seen that the equation of the charged condensed matter field gives the value of 
Tc while the equation of the scalar potential gives the behavior of order parameter at 
Tc. If the bulk theory has the same fundamental symmetry as the G-L theory, then 
the equation for the scalar potential leads to the mean field critical exponent. When 
the added cubic term breaks the symmetry, the critical exponent differs from 1/2. The 
so-called "non mean field" behavior in Ref. [22| is just the result of violation of the 
symmetry of the usual mean field theory like the G-L theory. In the present paper, 
we focus on the AdS^^/CFT'i situation, in which the dual superconductor is of 2 + 1 
dimension. The universal critical exponents of s-wave holographic superconductors in 
various spacetime dimensions have been studied in [^, in which the authors found that 
in the large N limit, the mean-field results are independent of dimension just like the 
G-L theory. It would be interesting to compute 1/iV-effects to see how the fluctuation 
will affect the "mean-field" behavior in the holographic superconductors!^. We obtain 
our conclusion by comparison of several holographic superconductors with mean field 
or non mean field critical exponents and the results are convincing. However, it is a 
better way to apply the holographic renormalization group ^ ^ analysis to these 
models to find the low energy effective theories of holographic superconductors. This 
may bridge the gap between the holographic models of superconductor and the G-L 
theory. 
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APPENDIX 



A. Variation method to solve the Sturm-Liouville problem 

To familiar the readers to the variation method to solve the Sturm-Liouville problem, 
we present some basic results of this method and all the derivations can be found in 
many text books such as . The Sturm-Liouville eigenvalue problem is to solve the 
equation 

ax ax 
with boundary condition 

k{x)y{x)y\x)\l = Q. (A.2) 
The Sturm-Liouville problem can be result to be a functional minimize problem: 

r . dx{k{x)y'{x)'^ + q{x)y{xY) 

dxp{x)y{xY 

where the minimal eigenvalue Aq and its eigenstate ^0(2^) can be obtained by variate the 
above equation. The n + 1 th eigensystem can also be obtained by variation eg ( |A.3|) 
with constrains below: 

h 

dxp{x)y*^{x)yi{x) = 5n,i (z = 0, 1, n - 1, ra). (A.4) 

then, the eigenvalue A„ satisfies: 

^ ^ Cdx{k{x)y\xf + q{x)y{xf) 
Iadxp{x)y{xy 

with complete eigenfunction {yn{x)}. 
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